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In this paper, we numerically investigate the length spectra and the low-lying eigenvalue spectra 
of the Laplace-Beltrami operator for a large number of small compact (closed) hyperbolic (CH) 3- 
manifolds. The first non-zero eigenvalues have been successfully computed using the periodic orbit 
sum method, which are compared with various geometric quantities such as volume, diameter and 
length of the shortest periodic geodesic of the manifolds. The deviation of low-lying eigenvalue 
spectra of manifolds converging to a cusped hyperbolic manifold from the asymptotic distribution 
has been measured by C,— function and spectral distance. 



I. INTRODUCTION 



Eigenmodes of the Laplace-Beltrami operator A on a Riemannian manifold (or orbifold) which carry information 
of both local geometry and global topology 

{J\ + E)uE^O (1) 

play a significant role in various kinds of physical systems. 

In the cosmological perturbation theory, one can interpret i?^/^ as the wavenumber k of an eigenfunction which 
characterises the scale of the metric or matter perturbation. In order to distinguish models with equivalent local 
geometry but non-equivalent global topology, it is crucial to study the behavior of low- lying eigenmodes. For spatially 
compact models, the lowest non-zero wavenumber ki (the first "excited" state) determines the maximum fiuctuation 
scale of the perturbation since fluctuations on scale less than are strongly suppressed. In compact flat 3-manifolds, 
is approximately equal to an inverse of the diameter(the maximum distance between two points) of the space. 
Therefore, if the space elongates in one dimension and shrinks in the remaining dimension keeping the volume constant, 
ki converges to zero. For compact hyperbolic(CII) 2-manifolds where one can deform the space continuously (without 
topology change) Rayleight's theorem implies that ki can be arbitrarily close to zero However, subtlety arises in 
the case of CH 3-manifolds. Because of the Mostow's rigidity theorem, one cannot deform the space continuously. 
Instead, we have a series of manifolds with different topology converging to a cusped manifold with finite volume 
(co-finite manifold). For CH 3-manifolds Af, orthonormal basis of L^(Af) may include supercurvarture modes fc< 1 
whereas there are not any such modes in the simply-connected hyperbolic 3-space H'^. If any supercurvature modes 
were present, the fluctuation property on large scales would be drastically altered. 

In quantum mechanics, ue can be interpreted as a wave function of a free particle at a stationary state with energy 
E. The statistical property of the energy eigenvalue E and the eigenfunction ue have been investigated for exploring 
the imprints of classical chaos in the corresponding quantum system(e.g. see [|] and other articles therein). Because 
any classical dynamical systems of a free particle in CH spaces are strongly chaotic (K-systems), the semiclassical 
behavior of the statistical property of eigenvalues and eigenfunctions in these spaces has been intensively studied [^-@ . 
It is the Gutzwiller trace formula jl^] that relates a set of periodic orbits(=geodesics) to a set of energy eigenstates 
and gives the semiclassical correspondence for classically chaotic systems. Interestingly, for CH spaces, the Gutzwiller 
trace formula gives an exact relation which had been known as the Selberg trace formula in mathematical literature 
[ p2| . The trace formula gives an alternative method to compute the eigenvalues in terms of periodic orbits. The poles 
of energy Green's function are generated as a result of interference of waves each one of which corresponds to a periodic 
orbit. Roughly speaking, periodic orbits with shorter length contribute to the deviation from the asymptotic eigenvalue 
distribution on larger energy scales. In fact, zero-length orbits produce Weyl's asymptotic formula. Because periodic 
orbits can be obtained algebraically, the periodic orbit sum method enables one to compute low-lying eigenvalues for 
a large sample of manifolds or orbifolds systematically if each fundamental group is known beforehand. The method 
has been used to obtain eigenvalues of the Laplace-Beltrami operator on CH 2-spaces and a non-arithmetic 3-orbifold 
P,|6UlO|]. However, it has not been applied to any CH 3-manifolds so far. 

To date, various numerical techniques have been applied to the eigenvalue problems for solving the Helmholtz 
equation (pi) with periodic boundary conditions (manifold case), Neumann and Dirichlet boundary conditions (orbifold 
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case). Eigenvalues of CH 2-spaces has been numerically obtained by many authors Eigenvalues of 

cusped arithmetic and cusped non-arithmetic 3-manifolds with finite volume have been obtained by Grunewald and 
Huntebrinker using a finite element method . Aurich and Marklof have computed eigenvalues of a non- arithmetic 
3-orbifold using the direct boundary element mcthod(DBEM) To date, computation of eigenvalue spectra is 

limited to a small number of 3-manifolds. The author has succeeded in computing eigenvalues of the Thurston 
manifold, the second smallest one using the DBEM and later the Weeks manifold, the smallest one in the known 
CH manifolds using the same method Cornish and Spergcl have also succeeded in calculating eigenvalues of 

these manifolds and 10 other CH manifolds based on the Trefftz method [T7| . 

In this paper, we study the length spectra and the low-lying eigenvalue spectra of the Laplace-Beltrami operator 
on a relatively large number of small CH manifolds which are obtained by a computer program "SnapPea" by Weeks 
p8[ . We analyse the fluctuating property of length spectra and check the accuracy of the first non-zero eigenvalues 
obtained by the periodic orbit sum method based on the trace formula. We also relate the low-lying eigenvalues to 
diameter, volume and length of the shortest periodic orbit and characterise the deviation of the spectrum from the 
asymptotic distribution. In Sec. II we briefly describe fundamental aspects of CH 3-manifolds which we will study. In 
Sec. HI we study the length spectra of CH 3-manifolds, especially we put an emphasis on their fluctuating property. 
In Sec. IV we derive an explicit form for computing the spectral staircase in terms of length spectra from the trace 
formula. In Sec. V we analyse the relation between the low- lying eigenvalues and several diffeomorphism-invariant 
geometric quantities, namely, volume, diameter and length of the shortest periodic orbit. In Sec. VI the deviation 
of the low-lying eigenvalue spectra from the asymptotic distribution for manifolds that have a region similar to the 
neighbourhood of a cusped point is measured by (— function and the spectral distance. 



II. HYPERBOLIC MANIFOLDS 



The discrete subgroup F of PSL{2,C) which is the orientation-preserving isometry group of the simply-connected 
hyperbolic 3-space H'^ is called the Kleinian group. Any CH 3-spaces (either manifold or orbifold) can be described 
as compact quotients M = W^/T. If we represent H'^ as an upper half space (xi, a;2, 0:3), the metric is written as 

where R is the curvature radius. In what follows, R is set to unity without loss of generality. If we represent a point 
p on the upper-half space, as a quaternion whose fourth component equals zero, then the actions of PSL[2, C) on 
U C U {00} take the form 

^ : p ^ p' — ^ ^ , , ad—bc=l, p = z + X3j, z = xi + X2I, (3) 
cp + d 

where a, b, c and d are complex numbers and 1, i and j are represented by matrices, 

'^(o-O' 

The action 7 is explicitly written as 

7 : U C U {00} ^ U C U {00}, 



, , , , I {az + h){cz + d) + acxi x^ \ 

where a bar denotes a complex conjugate. Elements of F for orientable CH manifolds are conjugate to 

/exp(V2 + *0/2) \ 

^\ exp(-V2 - i(/)/2) j 

which are called loxodromic if </) 7^ and hyperbolic if = 0. 

Topological construction of CH manifolds starts with a cusped manifold with finite volume Mc obtained by gluing 
ideal tetrahedra. Let us consider the case where Mc is topologically equivalent to the complement of a knot K or 
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FIG. 1. Volume spectra of CH manifolds computed by SnapPea. A prefix "m" in the labelling number represents a cusped 
manifold obtained by gluing five or fewer ideal tetrahedra. The numbers in the right side denote the volumes of the corresponding 
cusped manifold. 



link i(which consists of knots) in 3-sphere S'^ or some other closed 3-spaces. A surgery in which one removes the 
tubular neighborhood N oi K whose boundary is homeomorphic to a torus, and replace A'^ by a solid torus so that 
a meridian^ in the solid torus goes to {p, q) curvej^ on N is called {p, q) Dehn surgery. Except for a finite number of 
cases, Dehn surgeries on K always yield CH 3-manifolds which implies that most compact 3-manifolds are hyperbolic 
. SnapPea can perform Dehn surgeries which have made it possible to construct a large number of samples of CH 
3-manifolds. 

It is known that only a finite number of CH 3-manifolds with the same volume exist JTof . Hence the volume plays 
a crucial role in describing CH 3-manifolds. The key facts are: the volumes of CH 3-manifolds obtained by Dehn 
surgeries on a cusped manifold Mc are always less than the volume of Mc, CH 3-manifolds converge to Mc in the 
limit \p\, \q\ oo. As shown in figure |l|, the volume spectra are discrete but there arc many accumulation points 
which correspond to the volumes of cusped manifolds. The smallest cusped manifolds in the known manifolds have 
volume 2.0299 which are labeled as "m003" and "m004"in SnapPea. m003 and m004 are topologically equivalent to the 
complement of a certain knot in the lens space L^ ^ and the complement of a "figure eight knot" (figure respectively 
pO| . (3,-1) and (-2, 3) Dehn surgeries on m003 yield the smallest and the second smallest known manifolds, which are 
called the Weeks manifold (volume=0.9427) and the Thurston manifold (volume=0.9814), respectively. As \p\ and \q\ 
becomes large, the volumes converge to that of Mc- Similarly, one can do Dehn surgeries on m004 or other cusped 
manifolds to obtain a different series of CH manifolds. 



III. LENGTH SPECTRA 



Computation of periodic orbits (geodesies) are of crucial importance for the semiclassical quantization of classically 
chaotic systems which will be discussed in the next section. However, in general, solving a large number of periodic 
orbits often becomes an intractable problem since the number of periodic orbits grows exponentially with an increase 
in length. For CH manifolds periodic orbits can be calculated algebraically since each periodic orbit corresponds to 
a conjugacy class of hyperbolic or loxodromic elements of the discrete isometry group F. The conjugacy classes can 
be directly computed from generators which define the Dirichlet fundamental domain of the CH manifold 

Fet gi, {i = 1, N) be the generators and T be the identity. In general these generators are not independent. They 
obey a set of relations 



^ Given a set of generators a and b for the fundamental group of a torus, a closed curve which connects a point x in the torus 
with ax is called a meridian and another curve which connects a point x with bx is called a longitude. 
■^If C connects a point x with another point {pa + qb)x where p and q are co-prime integer, C is called a (p, q) curve. 
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FIG. 2. "Figure eight knot" (left) and the Dirichlet domain of a cusped manifold m004 viewed from two opposite directions 
in the Klein(projective) coordinates where geodesies and planes are mapped into their Euclidean counterparts (right). The 
two vertices on the left and right edges of the polyhedron which are identified by an element of the discrete isometry group 
correspond to a cusped point. Colors on the faces correspond to the identification maps. One obtaines the Dirichlet domain of 
m003 by interchanging colors on quadrilateral faces in the lower(or upper) right figure. 



which describe the fundamental group of M. Since all the elements of F can be represented by certain products of 
generators, an element g (zT can be written 

9^9ii9i2---9i„, (8) 

which may be called a ^^word^^ . Using relations, each word can be shorten to a word with minimum length. Fur- 
thermore, all cyclic permutations of a product of generators belonging to the same conjugacy class can be eliminated. 
Thus conjugacy classes of F can be computed by generating words with lowest possible length to some threshold 
length which are reduced by using either relations (R) among the generators or cyclic permutations of the product. 

In practice, we introduce a cutoff length Icut depending on the CPU power because the number of periodic orbits 
grows exponentially in I which is a direct consequence of the exponential proliferation of tiles (copies of the funda- 
mental domain) in tessellation. Although it is natural to expect a long length for a conjugacy class described in a 
word with many letters, there is no guarantee that all the the periodic orbits with length less than Icut is actually 
computed or not for a certain threshold of length of the word. 

Suppose that each word as a transformation that acts on the Dirichlet fundamental domain D. For example, D 
is transformed to D' = gD by an element g. We can consider gD's as tessellating tiles in the universal covering 
space. If the geodesic distance d between the center (basepoint) of D and that of gD is large, we can expect a long 
periodic orbit that corresponds to the conjugacy class of g. Tessellating tiles to sufficiently long distance d > Icut 
makes it possible to compute the complex primitive length spectra {Lj = Ij exp{i(f)j), \lj < Icut} where Ij is the real 
length of the periodic orbit of a conjugacy class with one winding number and is the phase of the corresponding 
transformation. We also compute multiplicity number m{lj) which counts the number of orbits having the same lj 
and 

In general, the lower limit of the distance d for computing a complete set of length spectrum for a fixed Icut is not 
known but the following fact has been proved by Hodgson and Weeks . In order to compute a length spectra of a 
CH 3- manifold (or 3-orbifold) with length less than /, it suffices to compute elements {g} satisfying 

d{x,gx) < 2cosh~^(coshi?coshZ/2), (9) 
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FIG. 3. The classical staircases N{1) for the Weeks manifold and the Thurston manifold with the asymptotic distribution 



where a: is a basepoint and R is the spine radiufjj. Note that there is a unique geodesic which lies on an invariant axis 
for each hyperbolic or loxodromic element. SnapPea can compute length spectra of CH spaces either by the "rigorous 
method" based on the inequality (H) or the "quick and dirty" method by setting the tessellating radius d by hand. 
I have used the former method for manifolds with small volume (< 1.42), but the latter method {d — Icut + 0.5) 
has been also used for some manifolds with large volume (> 1.42) since the tessellating radius given by the former 
method is sometimes so large that the computation time becomes too long. The detailed algorithm is summarized in 
appendix A. 

The asymptotic behavior of the classical staircase N{1) which counts the number of primitive periodic orbits with 
length equal to or less than I for CH 3-spaces can be written in terms of I and the topological entropy r ||2^ ] 

7V(0 - Ei(rO - l^oo. (10) 

Tt 

The topological entropy for Z3-dimensional CH spaces is given by r = Z? — 1. A larger topological entropy implies 
that the efficiency in computation of periodic orbit is much less for higher dimensional cases [10| . 

In figure^, the computed classical staircases {I cut — 7.0) are compared with the asymptotic formula for the smallest 
(Weeks) manifold and the second smallest (Thurston) manifold. For both cases, an asymptotic behavior is already 
observed at I ^ 3.5. 

Although the asymptotic behavior of the classical staircase N{1) does not depend on the topology of the manifold, 
the multiplicity number m{l) does. In fact, it was Aurich and Steiner who firstly noticed that the locally averaged 
multiplicity number 

<ni{l)>~ — ^ 9{h), (A^ = total number of terms) (11) 

l-Al/2<li<l+Al/2 

grows exponentially <m{l) >~ e'/^/Z as l^oo for arithmetic 2-spaces (manifolds and orbifolds) |2^j^,Q. since the 
length / of the periodic orbits are determined by algebraic integers in the form 2cosh(Z/2) — algebraic integei^. The 
failure of application of the random matrix theory to some CH spaces may be attributed to the arithmetic property. 
For non-arithmetic spaces, one expects that the multiplicities are determined by the symmetries (elements of the 
isometry group) of the space. However, in the case of a non-arithmetic 3-orbifold, it has been found that < m(l) > 



^Spine radius R is equal to the maximum over all the Dirichlet fundamental domain's edges of the minimum distance from 
the edge to the basepoint. Note that R is finite even for a cusped co-finite manifold. 

*For a two-dimensional space, the classical staircase has an asymptotic form N{1) ~ exp{l)/l. On the other hand, the classical 
staircase for distinct periodic orbits has a form N{1) ~ exp(//2) for arithmetic systems. Because m{l)dN = dN we have 
m{l) ~ exp(//2)// as I oo. 
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FIG. 4. Plots of locally averaged multiplicities (Al = 0.2) with one-sigma errors and the fitting curves for the Weeks manifold 
m003(-3,l) which is arithmetic (star) and a non-arithmetic manifold m004(l,2) (diamond). The fitting curves (aexp(Z)// for 
the former and exp(6Z)/(c/) for the latter) are obtained by the least square method using data 3.0 < I < 7.0. 



grows exponentially in the form e^^ /[cl) where b and c are fitting parameters This fact might implies that the 
symmetries of long periodic orbits are much larger than that of the space even in the case of non-arithmetic systems. 

For 3- manifolds, one can expect that the property of < m{l) > for 3-orbifold also holds. The locally averaged 
multiplicities (Al = 0.2) for the smallest twelve examples which include seven arithmetic and five non-arithmetic 
3-manifolds have been numerically computed using SnapPea. From figure ^, one can see that the difference in 
the behavior of < m{l) > between the non- arithmetic manifold and the arithmetic one is manifest. As observed in 
3-orbifolds, averaged multiplicities behave as 

<m{l)> =a — (arithmetic), (12) 



I 

exp bl 



(non-arithmetic), (13) 



where a depends on the discrete isometry group while b and c are fitting parameters. From table I, one observes that 
arithmetic manifolds having a larger symmetry group have a larger value of a. The growth rates for non-arithmetic 
manifolds (6 ~ 0.56) are always less than that for arithmetic manifolds (6 = 1) but nevertheless exponential. 



IV. PERIODIC ORBIT SUM METHOD AND SPECTRAL STAIRCASE 



Gutzwiller's periodic orbit theory provides a semiclassical quantization rule for classically chaotic systems. The 
theory is expressed in form of a semiclassical approximation {h 0) of the trace of the energy Green's operator 
(resolvent operator) Ge = (A + E)~^ in terms of the length of periodic orbits (geodesies) {Li] which is known as 
the Gutzwiller trace formula [ pH . For the dynamical system of a free massive particle on a CH space known as the 
Hadamard-Gutzwiller model, the periodic orbits give the exact eigenvalues, and the relation is no longer semi-classical 
approximation. In mathematical literature, the trace formula is known as the Selherg trace formula jl^ . In what 
follows we consider only orientable CH 3-manifolds (denoted as CH manifolds) (for general cases including orbifolds. 



see |10 ). The Selberg trace formula for a CH manifold M = /T (F is a discrete isometry group containing only 



hyperbolic or loxodromic elements) can be written as 

Tr{GE-GE-)^~^-^{p-p') 
47ri 
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TABLE I. Coefficients of the fitting curves which describe the average behavior of locally averaged multiplicities < m{l) > 
for arithmetic (A) and non-arithmetic (N) 3-manifolds. G denotes the isometry group (symmetry group). Fitting parameters 
a, b and c are obtained by the least square method using data 3.0 < I < 7.0. f For m007(3,l), the isometry group may be larger 
than D2. 



KdTo) I exp{-ipl{gr)) exp{-ip'l{gr)) 

4('nnsh l(r 



^ 4(cosh;(5^) - cos (l){gr)) ip 



ip' 



(14) 



where — E — 1, v{M) denotes the volume of M, ligr) is the (real) length of the periodic orbit of transformation 
gr G r. grg is a transformation that gives the shortest length of the periodic orbit ligro) which commutes with gr- 
(j){gT) is the phase of the transformation gr- The sum in (|l4| ) extends over V conjugacy classes 

{gr}:^{g'r\9'r = hgrh-\h(^T} (15) 

of hyperbolic ((/> = 0) or loxodromic elements (0 7^ 0). However the periodic orbit sum in ( [1^ ) which is known as 
MaaB-Selberg series converges at only complex energy such that Im p<—l and Im p' <—l. In order to obtain real 
eigenvalues, one needs to multiply the trace by some suitable analytic "smoothing" function h{q) that satisfies: 
{l):h{q) = h{-q)- 

{n):h{q) = 0{\q\-^^^) for (5 > as \q\ oo\ 

(iii):/i(q) is analytic in the strip |Im g| < 1 + e for e > 0. 

Multiplying ( p^ ) by q h{q)/{'K i) and integrating it over q from —00 to 00, one obtains the general Selherg trace 
formula, 

^f ,i.os.,M±.,M, 'VMl (16) 

where 

h{l) = ^ / dq h{q) cxp{-iql), (17) 

and h"{0) is the second derivative of h and pn denotes a wavenumber of the corresponding eigenmode. The sum in 
( p^ ) is absolutely convergent for any real eigenvalues i?„ = + 1 . One can obtain various functions of eigenvalues 
such as heat kernels and energy level densities from periodic orbits by choosing an appropriate "smoothing" function 

In order to obtain eigenvalues, a simple approach is to compute the spectral staircase 

00 

7v(i?) = ^^(£; -£;„), (18) 

n=0 



n=0 



7 



where Eq = 0, Ei, E2, ■ ■ ■ are the eigenvalues of the Laplacian and 9 is the Heaviside function To explore super- 
curvature modes M£;,0<£'<l,we choose the "smoothing" function of N{E) as 



hneip') = 77 ( 1 - Erf 



1 / . „ JE' - E 
2 



1/ Jp''^-p'^ 



which is real for E > 0. Note that hp{p') satisfies all the conditions (i) to (iii). By taking the limit e ^ 0, one obtains 
the spectral staircase 

oc 

N{E)^]imY,hj,M. (20) 

n=0 

Let us first estimate the behavior of the trace in ([l6|). From a straightforward calculation, the zero-length contri- 
bution can be written as 

KM^ --^exp(VA') 

^''(r(3/4))-iF(5/4, 1/2, //e^) + |4r(3/4)F(7/4, 3/2, //e^)), (21) 



where r(a;) is the Gamma function and F{a,b,z) is the confluent hypergeometric function. For a; =p*/e^ — > 00, it is 
asymptotically expanded as 



r{b-a 



^f^-'l'j: - ""^^ ~ (X)-" + 0(x-^)|, (22) 



k n=0 ) 

where [a)m = r(a + m)/r(a). From ([2l|) and (p2h, in the lowest order, we have the average part, 



^(p)-lim-^/.;:.(0) = fPb^|i, (23) 
- +0 27r " OTT^ 



that gives the dominant term in the Weyl asymptotic formula for p>> 1. 
Next, we estimate the oscillating term in ( [l^ ) 



In the long orbit-length limit I » 1 with p > 0, the integrand in Eq.(24) oscillates so rapidly that the dominant 
contribution comes from g~p, or —p where (g^ — p^)^/e^ can be approximately given as — 4(q — p^p^ /e^. Then ( |2^ ) 
can be written as 

Thus each periodic orbit corresponds to a wave with wavelength 2-k/I and an amplitude which is exponentially 
suppressed with an increase in Z or a decrease in p. For a finite subset of length spectra Ij < lent, the appropriate 
choice for the smoothing scale is given by e = ap^^^ for p^>0 where a depends on Icut since for a reasonable value of 
the proportional factor a all the contributions from periodic orbits with large length I > Icut can be negligible. For 
p^ <0, an optimal choice can be obtained from a numerical computation of Eq.(|2j) directly. Comparing the obtained 
smoothed spectral staircase with the one based on the computed "true" eigenvalues using the direct boundary element 
method (DBEM), it is numerically found that for ZcMt = 7.0, an appropriate smoothing scale is given by 



8 




02468 10 02468 10 

k k 

FIG. 5. The spectral staircases N{k) for the Weeks manifold and the Thurston manifold obtained by the DBEM are compared 

with the average parts N{k) (solid curves), namely, Weyl's asymptotic formula ( p^ ) and smoothed spectral staircase (dotted 

curves) obtained by the periodic orbit sum method(POSM) using all periodic orbits / < 7.0. 



manifold 


volume 


fci(DBEM) 


fci(POSM) 


m(fci) 


Afci/fci 


m003(-3,l) 


0.9427 


5.27 


5.10 
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0.03 


m003(-2,3) 


0.9814 


5.40 


5.34 
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0.01 


m007(3,l) 


1.0149 


5.29 


5.37 
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0.02 


m003(-4,3) 


1.2637 


4.58 


4.31 


2 


0.06 


m004(6,l) 


1.2845 


4.53 


4.35 


1 


0.04 


m004(l,2) 


1.3985 


4.03 


3.93 


1 


0.02 


m009(4,l) 


1.4141 


5.26 


4.84 


2 


0.08 



1/2 

TABLE II. The first (non-zero) wavenumbers fci = E-^ are calculated using the DBEM and using the POSM for seven 
smallest manifolds, fci's agree with relative accuracy Afei/fci = 0.01 — 0.08. The multiplicity number of the first eigenmode is 
calculated using the DBEM. 



/0.116fc2 + 0.184 fc + 1.2 (fc<l) , . 

^ ~ \ 0.832fci/2 + 0.668 (fc> 1) ' ^^^> 

where E — = + 1. The eigenvalues can be computed by searching E at which N{E) — 0.5 becomes positive 
integer. It should be emphasised that the precision of computation depends on the value of Zcut which determines the 
resolution scale in the eigenvalue spectra. 

In order to get eigenvalues from the smoothed spectral staircase, one must take into account the effect of the 
multiplicity number (degeneracy number) for each eigenvalue since the spectral staircase is smoothed on larger scales 
for degenerated modes. Therefore, the numerical accuracy becomes worse if the eigenmode has a large multiplicity 
number. Fortunately, the order of the symmetry group is not so large for a small manifold (volume< 3). For instance, 
of the twelve smallest examples, nine manifolds have a symmetry group with order 4. If one assumes that the 
multiplicity number is either 1 or 2 then the deviation Afc from a precise value is approximately given by 

e-ik + Ak)^\\ _ 1±0.5 ^27) 

For instance, if one uses a length spectrum / < 7.0 then ( p7| ) gives AA; = 0.30, 0.33 for fc = 5.0, 3.0, respectively. If one 
permits the multiplicity number as much as 6, then the expected precision becomes Afc = 0.49, 0.54 for fc = 5.0,3.0, 
respectively. 

We can see from table II that the first eigenvalues calculated by using length spectra I < 7.0 for some smallest known 
CH manifolds lie within several per-cent of those obtained by the DBEM. Note that the eigenvalues are also consistent 
with those obtained by the Trefftz method For two examples in which the first non-zero mode is degenerated, 
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the eigenvalues are much shifted to lower values owing to the smoothing effect. From figure |5|, one can see that the 
curves of the obtained smoothed spectral stairs cross the "true" stairs at almost half height. A slight deviation from 
the average part of the spectral staircase is caused by the interference of waves each one of which corresponds to a 
periodic orbit. 



V. FIRST EIGENVALUE AND GEOMETRICAL QUANTITIES 

The estimate of the first (non-zero) eigenvalue Ei = k\ of the Laplace-Beltrami operator plays a critical role in 
describing the global topology and geometry of manifolds. A number of estimates of Ei for n-dimensional compact 
Riemannian manifolds M using diffeomorphism-invariant quantities have been proved in mathematical literature. 

First of all, we consider the relation between the first eigenvalue Ei and the diameter d which is defined as the 
maximum of the minimum geodesic distance between two arbitrary points on M . Various analytic upper and lower 
bounds of Ei in terms of d have been known. Suppose M with Ricci curvature bounded below by ^L{L > 0). Cheng 
and Zhou proved that Ei satisfies 



El > max 



l7r^_ 1 
2d2 " 4 ' 



d2 



16 



^L,^exp{-CnVL^/2) 



(28) 



where C„ = max[V« — 1, V^] Another lower bound has been obtained by Lu Suppose that the Ricci 

curvature of M is bounded below as Rab ^ —Kgab, {K > 0) for a some real number K where gab is the metric tensors 
of M. Then Ei satisfies 



El > max 



d2 



8 K 8 
d^^Y'd^ 



exp 



d^K 



d2 



1 + -^Kin - 1) exp --^/Kin - 1) 



As for upper bounds, the following theorem has been proved by Cheng 
than (n — l)c, then we have 

Ei<Ei{Va{c,d/2)) 



(29) 



281. Suppose M with Ricci curvature larger 



(30) 



where Vn{c,r) denotes a geodesic ball with radius r in the n-dimensional simply-connected space with sectional 
curvature c and Ei is the first Dirichlet eigenvalue. Setting L = K — 2,n — 3, and c = — 1, we obtain the upper 



and lower bounds of Ei for CH 3-manifolds. For the upper bound, (30) gives a simple relation, d < where 
Vi — k\ — 1. The physical interpretation is clear: the wavelength Ai = 21: jvi of the lowest non-zero mode must be 
larger than the diameter. 

Now, we compare the first eigenvalues of 263 examples of CH 3-manifolds with volume less than 3 which have the 
length of the shortest periodic orbit Imin > 0.3 (the Hodgson- Weeks census |2^) and of 45 other examples obtained 
by Dehn surgeries {\p\ < 17, \q\ < 14) on a cusped manifold mOO30 with the analytic bounds. The diameter of a CH 
3- manifold is given by the supremum of the outradius^ over all the bascpoints, which has been numerically computed 
using the SnapPea kernel^. The numerical accuracy is typically Ad — 0.03 — 0.09 depending on the topology of the 
manifold. 

As shown in figure ^, the eigenvalues are well described by an empirical fitting formula Ai = Pd, or 



47r2 

^2^2- 



(31) 



Applying the least square method for the 263 manifolds in the Hodgson- Weeks census (group A), and 53 manifolds 
obtained by Dehn surgeries on m003 (group B), the best-fit values (3 = 1.70, 1.73 have been obtained for each group, 



^The Hodgson- Weeks census with volume less than 3 also includes 8 manifolds obtained by Dehn surgeries on m003. 
^The outradius at a basepoint x is equal to the minimum radius of the simply-connected ball which encloses the Dirichlet 
domain at x. 

would like to thank J. Weeks for providing me a code to compute the diameter using the SnapPea kernel. 
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1 1.5 ^ 2 2.5 

d 

FIG. 6. Diameter d versus ki for 263 examples of CH 3-manifolds with Imin > 0.3 and i; < 3 (group A, circle) and 52 examples 
that are obtained by performing Dehn surgeries (51 examples for \q\ < 10 plus 1 example (p, q) — (16, 13) group B, diamond) 
on a cusped manifold m003 with the best-fit curves for group B corresponding to ( |3l| ) (dashed curve) assuming fci(cusp) = 1 
and (^) (solid curve) assuming fci(cusp) = 0.1. The fitting curves also agree with the computed values for group A. The upper 
and lower thick curves denote the analytic bounds (ESl), (E9|) and (pQ). 



respectively. Note that /3 ^ 1.7 is slightly larger than the values 1.3-1.6 for 12 examples in the previous result by 
Cornish and Spergel |ll7|] . The deviation from the fitting formula (|3l] ) is found to be remarkably small (with one sigma 
error Aki = 0.18,0.19 for each group, respectively), which implies the existence of much sharper bounds. 



The empirical formula (31) asserts that no supercurvature modes {i.e. fci < 1) exist in the limit d ^ oo where the 
manifold converges to the original cusped manifold. However, cusped manifolds may have some supercurvature modes 
even for those with small volume. Instead of (p^), we consider a generalised empirical formula 



fci = y(fci(cusp))2-|-i^, (32) 

where fci(cusp) < 1 is the smallest wavenumber for the original cusped manifold. Although no supercurvature modes 
were observed in this analysis, the non-existence of such modes in the limit d oo was not confirmed since the 
numerical accuracy becomes worse for manifolds with small ki and large d. 

Next, we consider the relation between diameter d and volume v of the manifold. Since diameter is given by 
the supreme of the outradius (minimum radius of a sphere which circumscribes the Dirichlet domain) all over the 
basepoints, one expect that v is estimated as the volume of a sphere (in a hyperbolic space) with radius r somewhat 
smaller than d if there is no region which resembles the neibourhood of a cusp ( "thin part" ^) or equivalently Imin is 
sufficiently large. Suppose that r = ad with a < 1 then the volume of a sphere 

V = 7r(sinh(2ad) - 2ad) (33) 

gives the approximate value of a CH manifold with diameter d. The best fit value for a sample of 79 manifolds with 
Imin > 0.5 is a = 0.69. If a manifold has a "thin" part then the relation ( |3^ ) is no longer valid since the diameter 
becomes too long. Let Vc be the volume of a cusped manifold Mc (with only one cusp) and v{d) be the volume of a CH 
manifold M obtained by a Dehn surgery on Ale. In the limit d — > oo, one can show that the following approximation 
holds (see appendix B): 



A "thin" part is defined as a region where the injectivity radius(a half the minimum length of the periodic orbit) is short. 
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1 1.5 2 2.5 



FIG. 7. Diameter d versus volume v for 79 manifolds (Imt-n > 0.5) (star), 184 manifolds (0.5 > Imm > 0.3) (circle) and 53 
manifolds obtained by performing Dehn surgeries on m003 (group B, diamond) with the best-fit curves (^) (dashed-dotted 
curve) and (H) (dashed curve). The solid curve denotes the volume of m003. 



f T exp(-2(d- do)) \ 
v{d) = wc 1 J— , (34) 



1 

where 6 denotes a ratio of the volume of the complementary part to that of the "thin" part and dp is the diameter of 
the complementary part. For a sample of 41 manifolds with diameter longer than 1.65 in group B, the best fits are 
6 = 0.0 and do = 0.25. As shown in figure 0, the volume-diameter relation for CH manifolds with large Imin is well 
described by the fitting formula (p3|). As Imin becomes smaller, or equivalently, d becomes larger, a CH manifold M 
converges to the original cusped manifold Mc in which ( |3^ ) holds. 

Finally, we look into the relation between the first eigenvalue and the volume. For CH manifolds with sufficiently 
large Imin, @l and (||) give 

v{h)^7r{suM9{ki))~g{h)), g{k,) = -j==^^===. (35) 

PVikiY - (fci(cusp))^ 

To be consistent with the Weyl's asymptotic formula which is valid for ki » 1 

ki{v) = \/(^)Vl, (36) 

the fitting parameters should satisfy a/(3 — 3 ■ 4~^/^7r~^/'^ « 0.44 which well agree with the numerically computed 
values a/(3 — 0.69/(1.6 — 1.7) — 0.41 — 0.43 provided that fci(cusp) = 1. One can see from figure || that both (|35| ) 
and (^6|) give a good estimate of the first eigenvalue for globally "slightly anisotropic" manifolds with v < 3 and 
Imin > 0.5. For manifolds with large d, v and small Imin, (|3^) and (|3^ ) give incorrect estimates for fci. In that case, the 
effect of the curvature cannot be negligible. In contrast to compact flat spaces, the volume of a sphere in hyperbolic 
space increases exponentially as the radius increases. Assuming that the relation (|35| ) holds, for sufficiently globally 
"isotropic" CH manifolds (large Imin), ki is significantly larger than that for compact fiat spaces with the same volume 
even if one assumes that ki(cusp) ~ 0. However, for CH manifolds converging to the original cusped manifold Mc, 
the formula ( p5| ) has to be modified. If Imin is sufficiently small and d is large while keeping the volume finite then M 
has a "thin" part similar to the neibourhood of a cusp. Then one can use an asymptotic formula (B4h instead of (p3). 



k,{v) = J(fcr(cusp))2 + ^.(^^_i,(i_,/,^)/2)2 ^ (37) 
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— u . ^ — u .± u, u.± u.z u.o ± ± . z. 

Ak^ d 
FIG. 9. The length of the shortest periodic orbit Imin versus deviation Aki from the fitting formula 
is assumed(left) and Imin versus diameter d (right) for 263 CH manifolds (group A). 



where ki (cusp) = 1 
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FIG. 10. Spectral staircases for m003(9,7) (dashed-dotted staircase, left) , m003(16,13) (solid staircase, right) and 
m003(16,13) without the shortest periodic orbit (dotted staircase, right) are shown in comparison with the corresponding 
Weyl's asymptotic formula (solid curves) which have been numerically computed by the POSM using all periodic orbits with 
length I < 7. N{k) = k perfectly fits the staircase of m003(16,13) for fc < 5 (thick line, right). 



where Vc is the volume of Mc- As shown in figure ^, ki shifts to a smaller value for manifolds with smaller Imin 
which have larger d. In the limit M Mc, the length of periodic orbits of M also converges to that of Mc except 
for the shortest orbit whose length Imin goes to zero. From the general Selberg trace formula (p^), one can see 
that the wave which corresponds to the shortest orbit has a large amplitude ~ 1/ coshZmin with a long wavelength 
~ 2TT/lmin- Therefore, the presence of a very short periodic orbit results in the deviation of the energy spectrum 
at low-energy region (small k) from the asymptotic distribution. To confirm this, the spectral staircase using the 
length spectra of m003(16,13) but the shortest periodic orbit is removed has been computed. Note that m003(16,13) 
is very similar to the original cusped manifold m003 having d — 2.75 and Imin — 0.0086. As shown in figure |lO|, the 
computed spectrum staircase agrees well with Weyl's asymptotic formula. The computed spectrum may coincide with 
the one for a manifold in which the "thin" part is cut off. On the other hand, the spectral staircases for m003(9,7) 
and m003(16,13) which have small Imin deviate from the average part (=Weyrs asymptotic formula (^)) (figure |o|). 
These manifolds have a "thin" part which is virtually one-dimensional object. Let us remind that Weyl's asymptotic 
formula for a n-dimensional compact manifold M is given by 

U.nV{M)k'' ^ 

V ^ (27r)" ^ ^ 

where a;„ = 27r"/^/(nr(n/2)) is the volume of the unit disk in the Euclidean n-space and u(M) is the volume of 
M. If one assumes that the Weyl formula still holds for small k, then the spectrum in the low-energy region which 
corresponds to fluctuations on large scales for these manifolds can be approximately described by the Weyl formula 
for n = 1, N{k) oc k/ir. For m003(16,13), it is numerically found that N{k) = k provides a good fit for k < 5. In the 
next section, we will measure the deviation from the asymptotic distribution using the low-lying eigenvalue spectra. 



VI. SPECTRAL MEASUREMENT OF GLOBAL ANISOTROPY 



Among many possibilities, we should choose physically well-motivated quantities for measuring the global 
"anisotropy" in geometry in terms of eigenvalue spectra. First, we consider (^-function which is relevant to the cos- 
moki^ical microwave background anisotropy. The angular power spectra for CH universes are approximately written 

C^-E^. (39) 
i=i 



where Fi{k) can be approximated as a polynomial function of k. In order to measure the "anisotropy", we will define 
the following parameter, 

oo I oo 

1=1 ' i=l 

where kf are the eigenvalues of the Laplace-Beltrami operator on a CH 3-manifold M with volume v and fc^j are 
the eigenvalues obeying Weyl's asymptotic formula with volume v (^6|). Note that the zero- mode fco = is not 
included in the summation. Here we only consider the case s > 1 which ensures the convergence of the sum. The 
numerical result shows the clear difference between the 'slightly anisotropic" manifold m003(-3,l) (d=: 0.82, Imin = 0.58) 
, "somewhat anisotropic" manifold m003(9,7) (d = 2.11, = 0.028) and "very anisotropic" manifold m003(16,13) 
((i = 2.75, Zmin = 0.0086) (figure ^l|). The presence of the fluctuations on large scales in these "anisotropic" manifolds 
shifts the corresponding ^-function to a larger value. For the case in which the shortest periodic orbit is removed 
from the length spectra of m003(16,13), one can see that the spectrum coincides with that obeying Weyl's asymptotic 
formula. 
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FIG. 11. Spectrum measurements of the geometry. The figure in the left shows A(s)=("(s)/("i„(s) for 3 examples of CH 
3-manifolds ,m003(9,7) (dashed-dotted curve), m003(16,13) (solid curve), m003(-3,l) (dashed curve) and m003(16,13) where 
the shortest periodic orbit is removed (dotted curve). The figure in the right shows ds{n) = ^ X^^Li In \ {ki/kn,i + k^i/ki) for 
the same examples where n is the cutoff number in the summation. 



Next, we consider the spectral distance ds proposed by Seriu which measures the degree of semi-classical quantum 
decoherence between two universes having one massless scalar field pO| , 

..|M,Ml.l|,„i(|4), (41) 

where fc? and kf are the eigenvalues of the Laplace-Beltrami operator on a compact n-manifold M and M, respectively. 
Here we choose eigenvalues fc^^ as kf. In practice we introduce a cutoff in the summation. It is numerically found 
that the summation converges rapidly for the 3 examples, namely, m003(-3,l), m003(9,7) and m003(16,13). The 
contribution of the first several terms dominates the summation (figure ^) . The result implies that a universe having 
a spatial geometry m003(9,7) or m003(16,13) semiclassically decoheres with a universe having a spatial geometry 
m003(-3,l) (figure |l|). 
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FIG. 12. Plots of a Dirichlet domain of the Weeks manifold m003(-3,l) (left) and that of m003(16,13) (right) viewed from 
two opposite directions in the Klein coordinates. The Dirichlet domain of m003(16,13) is quite similar to that of the original 
cusped manifold m003 (see figure m003(16,13) has a "thin" part which is similar to the neibourhood of a cusp. The colors 
on the faces correspond to the identification maps. 



VII. SUMMARY 



In this paper, the length spectra and low-lying eigenvalue spectra of the Laplace-Beltrami operator on small CH 
3-manifolds have been numerically investigated. CH 3-manifolds are relevant to a various kinds of physical systems. 
For instance, computing the CMB anisotropy in CH models is one of the key issue which have been investigated 
for several models [^,^,^ 33 1. In order to fully understand to what extent CH models are constrained by the 
observations, it is necessary to explore the property of low-lying eigenvalues and eigenmodes for a large number of 
manifolds (or orbifolds). 

First, The length spectra for a total of 308 CH 3-manifolds (volume less than 3) have been successfully computed 
using the SnapPea kernel up to the length / = 7.0. The asymptotic behavior in the classical staircase is found to be 
consistent with the known analytical formula which does not depend on the topology or symmetry of the manifold. 
Regarding the symmetry of the length spectra, it is well known that the arithmetic structure breaks a generic feature: 
the locally averaged multiplicity number grows exponentially as I is increased. However, the exponential behavior in 
the multiplicity number which has been noticed for a non-arithmetic 3-orbifold has been confirmed for the smallest 
five examples of non-arithmetic 3-manifolds, though the rate is smaller than that for the arithmetic manifolds. This 
may be related to the "hidden symmetry" which is the symmetry of finite sheeted covers of the manifold (which can 
be tessellated by the copies of the fundamental domain of the manifold) but not of the manifold itself. 

Next, the trace formula has been applied to these 3-manifolds for computing low-lying eigenvalues using the length 
spectra (POSM). Consistency with those obtained by the DBEM has been confirmed for several manifolds. It is 
numerically found that these manifolds do not have any supercurvature modes. It seems that the manifolds which 
supports supercurvature modes are either having large volume or large diameter. In order to confirm this, further 
investigation of eigenvalues for cusped manifolds is necessary. 

Thirdly, the first eigenvalues are compared to a various diffeomorphism-invariant quantities, namely, diameter, 
volume and the shortest length of the periodic orbits. The numerical results imply the existence of much shaper 
bounds for the first eigenvalues in terms of diameter. Some fitting formulae have been introduced and their validity 
has been checked. We have seen that CH 3-manifolds can be roughly divided into two categories: "slightly anisotropic" 
and "almost anisotropic" ones. The former has not any very short periodic orbits while the latter has. For example, 
manifolds which are very similar to the original cusped manifold are belonging to the latter category. It is found that 
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the deviation of the spectrum from the Weyl asymptotic formula for these manifolds is conspicuous even for manifolds 
with small volume. 

Finally, the global "anisotropy" in the spatial geometry has been measured by (^-function and the spectral distance 
for 3 examples of CH 3-manifolds. These measurements give a clear indication of the presence of a "thin" part in the 
manifold in terms of eigenvalue spectra. In other words, the physical quantities (the angular power spectra in the 
CMB or decoherence between two universes) are greatly affected by the globally anisotropic structure in the spatial 
geometry. 
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APPENDIX A: 



For a given Dirichlet domain D, SnapPea computes 

1. Neighboring copies of D (tiles) recursively and stores the corresponding elements g of the discrete isometry 
group r. If it has already been found out, it is discarded. The computation proceeds until for all the neighborhoods 

of gD, d{x,hx) > cosh""'^ (cosh i? cosh ^/2) satisfies where hD's are neighborhoods of gD{R is the spine radius and h 
is an element of F). Since there is no g other than identity where all h's satisfy d{x,hx) > d{x,gx), this algorithm 
will not miss any tiles gD where d{x,gx) < 2cosh~^(coshi?coshZ/2). 

2. A list of geodesies for all g's = {g} where l:the real part of length is not zero and less than I; 2:the distance 
from x to the geodesic is at most R. 

3. The conjugacy class g' = hgh~^ or its inverse for each g where h is an element of {g}. If an identical com- 
plex length is found, the complex length which corresponds to g' is omitted from the list of geodesies. If the complex 
length of (/ is conjugate to that of g~^, the geodesic is topologically a mirrored interval, otherwise it is a circle. 

4. Multiplicity of geodesies. If a pair of geodesies with two complex lengths being identical within an error range is 
found, the multiplicity number is increased by one. 



APPENDIX B: 



Suppose a CH manifold M which resembles the original cusped manifold with one cusp. Let us divide 
into two parts, the neighbourhood of a cusped point Kc, and the complementary part Kco- Similarly, one can divide 
M into K and Kq where Kco ~ Kq and K corresponds to a "thin" part. Since the neighbourhood of a cusp is 
represented as a "chimney" (but having infinite length) in the upper half space coordinates {xi,X2,Xs), K can be well 
approximated by an elongated box defined by {—Ax/2 < xi < Aa;/2, — Aa;/2 < X2 < Ax/2,xsq < xs < xs\). Then 
the physical lcngth(=diameter) doiK m the direction xz is given by d = ln(a;3i/a;3o). On the other hand, the volume 
V oi K satisfies 



(Ax)2 / 1 1 



(X30)^ (2^31)^ 



(Bl) 



which gives the ratio of the volume v to the volume of Voo of Kc, v/voo = 1 — (2^30 )^/ (2:31)^ = 1 — cxp(— 2rf). If we 
approximate the diameter d oi M as d = d + do where do is the diameter of Ko, then we finally have the ratio of the 
volume V = vo+v of M to the volume Vc = vo + Voo of M^, 

V exp(-2(rf - rfo)) 

— = 1 . ' -, d = vo v^ (B2) 

Vc + 1 

where vq denotes the volume of Kq. 
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